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On . Abstract 



The interaction of hydrogen with many transition metal surfaces is characterized by a 

D ' coexistence of activated with non-activated paths to adsorption with a broad distribution 

Mh , of barrier heights. By performing six-dimensional quantum dynamical calculations using a 

potential energy surface derived from ab initio calculations for the system H2/Pd(100) we 

show that these features of the potential energy surface lead to strong steering effects in the 

^ . dissociative adsorption and associative desorption dynamics. In particular, we focus on the 

r^ I coupling of the translational, rotational and vibrational degrees of freedom of the hydrogen 

^— V ' molecule in the reaction dynamics. 

(N 

O, 

g;;; ■ 1 Introduction 

^ ; It is a long-term goal in surface science to understand catalytic reactions occuring at surfaces. 

Obviously, the single steps of these often rather complicated processes are more effectively studied 
'-^ ■ at simple systems. In particular, the dissociative adsorption and associative desorption of hydrogen 
^ , on metal surfaces has served as a benchmark system, both experimentally and theoretically (see, 
Q ', e.g., Refs. |l|, ^, §] and references therein). Since the mass mismatch between hydrogen and 
j> ' a metal substrate is rather large, the crucial process in the dissociative adsorption for these 
kp particular systems is the conversion of translational and internal energy of the hydrogen molecule 
j!^ ■ into translational and vibrational energy of the adsorbed hydrogen atoms. If in addition no surface 

rearrangement occurs upon adsorption, the substrate degrees of freedom can be neglected and the 

dissociation dynamics can be described in terms of potential energy surfaces (PES) which take 

only the molecular degrees of freedom into account. 

The PES for the dissociative adsorption of a diatomic molecule neglecting the substrate degrees of 
freedom is still six-dimensional. These PESs now become available by elaborate density-functional 
calculations @, H, |^, |^, ||. However, in order to understand the reaction dynamics one has to per- 
form dynamical calculations on these potentials. Because of its light mass hydrogen has to be 
described quantum mechanically. Only recently it has become possible to perform dynamical cal- 
culations of the dissociative adsorption and associative desorption where all degrees of freedom of 
the hydrogen molecule are treated quantum mechanically |^. These calculations have established 
the importance of high-dimensional effects in the reaction dynamics. 

For example, molecular beam experiments of the dissociative adsorption of H2 on various transi- 
tion metal surfaces like Pd(lOO) ^, Pd(lll) and Pd(llO) fTT]], W(lll) 0, W(IOO) |T|, 0, [P 



W(100)-c(2x2)Cu |jT5[ and Pt(lOO) [|lg] revealed that the sticking probability in these systems 
initially decreases with increasing kinetic energy of the beam. High- dimensional quantum dynam- 
ical calculations have shown that steering effects can cause such an initial decrease in the sticking 
probability |P, |T^ ; it is not necessarily due to a precursor mechanism, as was widely believed. In 
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Figure 1: Contour plot of the PES 
along a two-dimensional cut through 
the six-dimensional coordinate space of 
H2/Pd(100). The inset shows the ori- 
entation of the molecular axis and the 
lateral H2 center-of-mass coordinates. 
The coordinates in the figure are the H2 
center-of-mass distance from the sur- 
face Z and the H-H interatomic dis- 
tance d. The dashed line is the opti- 
mum reaction path. Energies are in eV 
per H2 molecule. The contour spacing 
is 0.1 eV. 



0.5 1.0 1.5 2.0 2.5 
H-H distance (A) 

this paper we briefly review the most important aspects of the steering mechanism, but mainly 
focus on the combined influence of the translational, rotational and vibrational degrees of free- 
dom of the hydrogen molecule on the adsorption and desorption process. In the next section the 
theoretical background will be introduced before the results of the dynamical calculations will be 
discussed. 



2 Theoretical background 

The potential energy surface of H2/Pd(100) has been determined using density-functional theory 
together with the generalized gradient approximation (GGA) [|l^ and the full-potential linear 
augmented plane wave method 



rU| , ^. Ab initio total energies have been evaluated for more 
than 250 configurations [^] and have been parametrized in a suitable form for the dynamical 
calculations |^. 

Figure |1] shows a cut through the PES of H2/Pd(100), where the most favourable path towards 
dissociative adsorption is marked by the dashed line. For this path there is no energy barrier 
hindering dissociation, i.e., the adsorption is non-activated. However, the majority of pathways 
towards dissociative adsorption has in fact energy barriers with a rather broad distribution of 
heights and positions, as the detailed total-energy calculations showed Q], i.e. the PES is strongly 
anisotropic and corrugated. This has important consequences, as will be shown below. 
The quantum dynamics is determined in a coupled-channel scheme within the concept of the 
local reflection matrix (LORE) [^. This numerically very stable method is closely related to the 
logarithmic derivative of the solution matrix and thus avoids exponentially increasing outgoing 
waves which could cause numerical instabilities. The reported calculations, which take all degrees 
of freedom of the hydrogen molecule into account, are still only possible if all symmetries of the 
scattering problem are utilized. 
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Figure 2: Sticking probability versus kinetic energy for a H2 beam under normal incidence on a 
Pd(lOO) surface. Theory: orient at ionally averaged sticking probability for different initial rota- 
tional quantum numbers < jj < 4 of the incoming molecular beam. The molecular beams are 
assumed to have an energy spread of AE/ Ei = 2Av/vi = 0.2 [|I^ {Ei and Vi are the initial kinetic 
energy and velocity, respectively). Experimental data (circles) are from ref. 
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3 Results 

Figure ^ presents six-dimensional quantum dynamical calculations of the sticking probability as 
a function of the kinetic energy of a H2 beam under normal incidence on a Pd(lOO) surface for 
different initial rotational states averaged over the azimuthal quantum numbers. In addition, the 
results of a molecular beam experiment are shown [|T0|. Quantum mechanically determined sticking 
probabilities for hydrogen at surfaces with an attractive well exhibit an oscillatory structure as a 
function of the incident energy [^, |I^, ^, ^ , reflecting the opening of new scattering channels and 
resonances [^, ^. These structures are known for a long time in He and H2 scattering [|^ and 
also in LEED P3|. For H2/Pd(100), however, measuring these oscillations is a very demanding 



task. They are very sensitive to surface imperfections like adatoms or steps |2^. Since we do not 
focus on these oscillations here, for the results of Fig. § we have assumed a velocity spread of the 
incoming beam typical for the experiment |T^ so that the oscillations are smoothed out. 
The initial decrease of the sticking probability found in the experiment is well-produced for 
molecules initially in the rotational ground state ji = 0. The high sticking probability at low 
kinetic energies is caused by the steering effect: Slow molecules can very efficiently be steered to 
non-activated pathways towards dissociative adsorption by the attractive forces of the potential. 
This mechanism becomes less effective at higher kinetic energies where the molecules are too fast 
to be focused into favourable configurations towards dissociative adsorption. This causes the ini- 
tial decrease of the sticking probability. If the kinetic energy is further increased, the molecules 
will eventually have enough energy to directly traverse the barrier region leading to the final rise 
in the sticking probability. 
The dynamical origin of the steering effect is refiected in the dependence of the sticking probability 



0.20 

0.18 

> 0.16 

^ 0.14 

I 0.12 

0.10 

"■1— « 

1 0.08 
c 0.06 
^ 0.04 

0.02 



0.00 



X 


equilibrium y'^ 


- rotationally averaged ,-''' j/^' 


- Jf = ° X' ^/^"": 




- j! = 4 / ^^^^^^^ - 


'-' -^x^^^" 


^:y^ -^ 




-''' .cy^-^^" 


-' -^;^^" 


'i^<<^-' 


,-' yji!^ 




^^ •iSp**^'^ 


-■^^^^^^ - 


1.1,1,1.1. 



200 300 400 500 600 700 

Surface temperature (K) 



800 



Figure 3: Mean kinetic energy of hydrogen molecules desorbing perpendicurlarly from a Pd(lOO) 
surface for different final rotational states as a function of the surface temperature. Thin dotted 
line: result for molecules in equilibrium with the surface temperature; thick solid line: summed 
over all final rotational states; dashed and dash-dotted lines: results for final rotational states 
Jf = 0, 2, 4. 



on the initial rotational state of the impinging molecules. Rapidly rotating molecules will rotate 
out of favourable configurations towards dissociative adsorption, i.e. the steering mechanism will 
be suppressed. This rotational hindering is demonstrated in Fig. §: at low kinetic energies the 
sticking probability decreases with increasing initial rotational quantum number ji. At higher ki- 
netic energies, where direct dissociative adsorption becomes dominant, rotational hindering is less 
important. As a consequence in this energy range the sticking probability is almost independent 
of the initial rotational state. 

The predicted rotational hindering in the sticking probability at low kinetic energies |^ has in 
the meantime been confirmed by experiment |^. The influence of the steering mechanism on 
the dissociation dynamics can also be demonstrated at the time-reverse process, namely the asso- 
ciative desorption. Invoking the principle of microscopic reversibility, any degree of freedom that 
hinders dissociative adsorption will in thermal associative desorption be less populated than ex- 
pected for thermal equilibrium with the surface temperature Tg. At low kinetic energies increasing 
the initial rotational state as well as increasing the kinetic energy hinders dissociative adsorption. 
In desorption this leads to the so-called rotational cooling, which was found experimentally |2^ 
and confirmed theoretically 0, as well as to translational cooling. The translational cooling is 
demonstrated in Fig. |^ where the mean kinetic energy in desorption as a function of the surface 
temperature for different final rotational states is shown. The thin dotted line corresponds to 
thermal equilibrium with the surface temperature. At low surface temperatures the mean kinetic 
energy of molecules desorbing in rotational states jf < 2 is indeed smaller than the equilibrium 
value. Translational cooling is also found if the mean kinetic energy is summed over all rotational 
states at low surface temperatures since at these temperatures the desorption flux is mostly pop- 
ulated by molecules in low rotational states. For large rotational quantum numbers the steering 
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Figure 4: Sticking probability versus kinetic energy for a H2 beam under normal incidence on 
a Pd(lOO) surface. Orientationally averaged sticking probability for different initial rotational 
quantum numbers < jj < 4 for molecules initially in the first excited vibrationally state Vi = 1. 
The molecular beams are assumed to have an energy spread of AE / Ei = 2Av/vi = 0.2. 



effect in adsorption is suppressed (see Fig. |^), consequently for jj = 4 only translational heating is 
found in desorption. Equivalently, at higher kinetic energies the sticking probability rises regard- 
less of the rotational state which leads to translational heating at higher surface temperatures for 
all final rotational states. 

However, Fig. |^ shows that the mean kinetic energy in desorption does not differ significantly 
from the thermal equilibrium value. State-resolved measurements of the mean kinetic energy in 
desorption were done in the surface temperature range 450 K <Tg < 850 K [Q. For H2 desorption 
the experiments do not show any translational heating, within the error bars the experiments are, 
however, consistent with our results. For D2 desorption on the other hand, translational heating 
has been found |^D[ in good agreement with our H2 calculations. 

Interestingly, strong vibrational heating in D2 desorption from Pd(lOO) has been found experi- 
mentally 1^. Usually this is associated with a so-called late barrier to adsorption as found in the 
system H2/CU 0, |^, |]. Since the adsorption of H2/Pd(100) is apparently non-activated p, |10|, |30| , 
the origin of the vibrational heating has been controversely debated [p2i |33| , ^ [35[| . Just recently 



we have shown that the vibrational effects in the system are caused by the lowering of the H-H 
vibrational frequency during the dissociation and the multi-dimensionality of the potential energy 
surface p6[ . 

From Fig. |l] it is apparent that the H-H vibrational frequency is strongly lowered during the dis- 
sociation. This is also true for other pathways to dissociative adsorption. A detailed analysis of 
the dynamics has further shown [Q that the H-H vibrations follow the change of the vibrational 
frequency almost adiabatically during the dissociation. This leads to an effective vibrational- 
translational energy transfer, or equivalently, to an effective lowering of the PES while the molecules 
approach the surface. The crucial point is that this lowering depends on the vibrational state: the 
higher the vibrational quantum number is, the larger the effective lowering is. Thus molecules in 



the first excited vibrational state experience a lower barrier range than molecules in the vibrational 
ground state. 

In fact, for molecules in the first excited vibrational state the PES due to the lowering is so attrac- 
tive that the steering mechanism is hardly operative. This is demonstrated in Fig. ^ where the 
sticking probability for molecules initially in the first excited vibrational state I'i = 1 for different 
initial rotational states j, is shown. For ji = only at very low kinetic energies Ei < 0.02 eV there 
is a significant decrease in the sticking probability with increasing kinetic energy. At Ei = 0.1 eV 
the sticking probability for vibrationally excited molecules is about 2.5 times larger than for 
molecules in the vibrational ground state (compare Fig. ^). In desorption this leads to a popu- 
lation of the first excited vibrational state which is two times larger than expected for thermal 
equilibrium with the surface temperature at Tg = 700 K p6| . 

The fact that for vibrationally excited molecules already at small kinetic energies direct disso- 
ciative adsorption is the dominant mechanism is also reflected by the dependence of the sticking 
probability on the initial rotational state for vibrationally excited molecules. Fig. ^ shows that 
this dependence is much less pronounced than for molecules in the vibrational ground state. Due 
to the effectively lowered potential vibrationally excited molecules are so fast that the rotational 
hindering does not play any crucial role. 

So far all the results presented were averaged over the azimuthal quantum number m, i.e, the 
results corresponded to orient at ionally averaged properties. Now the H2/Pd(100) PES is strongly 
anisotropic with regard to the molecular orientation. The most favourable configuration towards 
dissociative adsorption is with the molecular axis parallel to the surface. Molecules that hit the 
surface in an upright position cannot dissociate, they are reflected back into the gas-phase. It 
is true that quantum mechanics does not allow for non-rotating, oriented molecules in the gas- 
phase, however, rotating molecules can show a preferential orientation. Molecules with azimuthal 
quantum number m = j have their axis preferentially oriented parallel to the surface. It has been 
shown theoretically that these molecules rotating in the so-called helicopter fashion dissociate 
more easily than molecules rotating in the cartwheel fashion {m = 0) with their rotational axis 
preferentially parallel to the surface 0. 

Experimentally it is hard to align a molecular beam of hydrogen. Again one can study the time- 
reverse process, the associative desorption. By laser-induced fluorescence (LIE) it is possible to 
measure the rotational alignment parameter Aq (j) ^^, which is given by 

Aq '{j) corresponds to the quadrupole moment of the orientational distribution and assumes values 
of— 1 < Aq (j) < 2. Molecules rotating preferentially in the cartwheel fashion have an alignment 
parameter Aq (j) < 0, for molecules rotating preferentially in the helicopter fashion Aq (j) > 0. 
Figure ^ shows a comparison between experiment [^ and theory of the rotational alignment 
of hydrogen desorbing from a Pd(lOO) surface at a surface temperature of Tg = 690 K. Note 
that the experiments were done for D2, while the calculations were performed for H2, hence the 
comparison has to be done with caution. Still the agreement for j < 6 is quite satisfactory. Indeed 
the molecules desorb preferentially with their molecular axis parallel to the surface thus reflecting 
the anisotropy of the PES. Due to computational restrictions the rotational alignment could only 
be calculated for j < 6. For j = 7 and j = 8 the experiments show a vanishing alignment 
within the error bars which is rather surprising. It was argued that for these high rotational state 
an energy transfer from the rotations to the reaction coordinate could suppress the effect of the 
potential anisotropy pH|. Certainly these results deserve further clariflcation. 
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Figure 5: Rotationally alignment of hydrogen molecules desorbing from a Pd(lOO) surface. Boxes: 
experimental results for D2 p8[. Circles: 6-D calculations for H2. 



4 Conclusions 

In conclusion, we reported a six-dimensional quantum dynamical study of dissociative adsorption 
and associative desorption of H2/Pd(100). We have briefly reviewed the steering mechanism, in 
particular in the rotational coordinates, which leads to the initial decrease of the sticking prob- 
ability with increasing kinetic energy. Due to the multi-dimensionality of the relevant potential 
energy surface the translational, rotational and vibrational degrees of freedom of the hydrogen 
molecule are strongly coupled during the dissociation. Our results establish the importance of a 
high- dimensional dynamical treatment in order to understand reactions at surfaces. 
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